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Abstract 

The goal of the present paper is to present a new approach to 
the construction of asymptotic (approximating) solutions to parabolic 
PDE by using the characteristics. 



1 Forward in time motion 
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The goal of the present paper is to present a new approach to the construc- 
q ■ tion of asymptotic (approximating) solutions to parabolic PDE by using 

the characteristics. This approach allows one to construct global in time 
solutions not only for the usual Cauchy problems but also for the inverse 
problems. We will work with Kolmogorov-Feller-type equations with diffu- 
sion, potential, and jump terms. The equation under study has the form: 

l 

du _ / 2 d \ 

u(x,t,e)\ t=0 = e- So ^ £ <p (x), 

where P(x, £) is the symbol of the Kolmogorov-Feller operator, e — > +0 is a 
small parameter characterizing the frequency and the amplitude of jumps of 
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the Markov stochastic process with transition probability given by P(x,£). 
To be more precise, we bear in the mind the following form of P(x,£): 



where A(x) is a positive smooth matrix, p,{x, dv) is a family of positive 
bounded measures smooth with respect to x such that 



and B and V are smooth in x (more precise conditions see below). The 
construction of forward in time global asymptotic solution to equations of 
this type was developed by V. Maslov, [U QjJJ [23], for a version of this con- 
struction, see also in [201 E21 [21] . Maslov's approach is based on ideas similar 
to those used in his famous canonical operator construction (or in Fourier 
integral operators theory). This construction is based on some integral rep- 
resentation and is not suitable for constructing backward in time solutions. 
Another approach to the global asymptotic solution construction was sug- 
gested in [2] and is based on the construction of generalized solutions to 
continuity equation in a discontinuous velocity field [1|. 

We assume that the class of solutions under study admits the following 
limits: 

(1) - logarithmic pointwise limit lim £ __ > o( — £lnu). We denote this limit by 
S = S(x,t) and assume that it is a piecewise smooth function with bounded 
first-order derivatives and a singular support in the form of a stratified 
manifold M. 

(2) - weak limit of the expression exp (2S / e)u 2 . We denoted it by p 
and assume that p is the sum of the function (p r eg) smooth outside M and 
the Dirac 5-function on M. Note that here we deal with the limit in the 
weighted weak sense! 

If S and p are smooth function, then the following representation is true 
(in the usual sense) 



P(x, = (A(x)£, £) + V(x) + (B(x),0 + 





i = 1, ... ,n 



u = exp (-5/e)v^ reg (l + o{e)). 



Example 1: 



u\ t =o = exp (-So/e)ip (x) 
where Sq > is a smooth function, ip° £ C^°. 
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Here a WKB-like approach can be used (Yu. Kifer, [3j; V. Maslov, [4, 
It gives an asymptotic (approximating) solution in the form (cf. [7]) 

Uas = Uas(x,t,£) = exp (S(x,t)/s)((p (x,t) H h e k (f k (x,t)) 

for arbitrary k. Here S(x, t) is the solution to the Cauchy problem for the 
Hamilton-Jacobi equation 

S t + P(x,VS)=0,S t=0 = S (x), (1.2) 

and (po = (fo(x,t) is the solution to the transport equation 

<Pot + (V s P(x, V5), Vv9 ) + ytr(P K (x, V5)5 XiE = 0, (1.3) 

St=o = S (x). 

Both of the solutions S and y? are defined via solutions of the Hamilton 
system 

x = V^P(x,p), x\ t =o = a,p = -V x P(x,p), p\ t =o = VS (a). (1.4) 
They are smooth while 

Da r 

There are symplectic geometry objects corresponding to this construc- 
tion: 

(1) the phase space R 2 x n p = M™ x R"; 

(2) the Lagrangian manifold G N* 1 , 

A° = (x = a,p = VS (a)), 

where g l p is a shift mapping along the Hamiltonian system trajectories; 

(3) the projection mapping ir : y 1R" with Jacobi matrix t^j. 



The main assumption that is required is the following one. 

The trajectories of the Hamilton system form a manifold of the phase 
space (at least in the area of the phase space under study). 

Let Dx/Da ^ for t 6 [0,T], then we have the following statement 
(V. Maslov, mil]; V. Danilov, [201 122J EE] ) : 
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Theorem 1. The inequality 

ReP(x,p + irj) < P(x, p), n G R n 

is necessary and sufficient for the following estimation to hold: 

\\exp(S(x,t)/e)(u as - u)\\ C (Rn x ) < C M £ M , 

where M = M(k) — > oo as k — > oo and the function S is a solution to the 
Hamilton- J acobi equation with Hamiltonian P(x,p). 

It is easy to verify that the function P = P(x, £) - the symbol introduced 
above - satisfies the inequality mentioned in the theorem. 
Example 2: 

-eu t + e 2 u xx = 0, u\ t =o = exp (-So/e)f . 
The corresponding system is: 

St + (S x ) 2 = 0, (Hamilton Jacobi equation), 



(fot + 2S x ipo x + S xx (fo = 0, (transport equation) 



x = 2p, x\ t =o = a, 
dS 

p = 0, p\t=o = -q — i (Hamilton system). 



Its solution has the form 

dS 

x = a + 2t— — , p = p\t=o 
oa 

and 

Dx/Da = l + 2t—4. 

(i) 

(ii) 

For t > t* in case (ii), we get A} of the shape plotted in Fig. 3. 
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Figure 1: Solution a = a(x,t) exists for all t 




Figure 2: Solution a = a(x,t) exists for t < t* = max | 2S'q qq , | 1 




Figure 3: Lagrangian curve for t > t*. 

One can see that in the case under study there are three values of S at 
the point x. This means that we can present an asymptotic solution near 
this point in the form of a linear combination: 



u 



E 



CjUj, 



(1.5) 
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where each of the functions Uj = exp(—Sj/e)ipj, j = 1,2,3, satisfies the 
equation with the same accuracy. But the functions themselves are not 
equivalent in contrast to the hyperbolic case. 
For example, it is clear that if the inequality 

Si(x,t) > S 2 (x,t) 

holds at a certain point x, then the "WKB" solutions U\ and u 2 at the point 
x satisfy the relation 

= e -5 1 (^)/ £(/92(S)t)(e -(5 1 -S 2 )/ ViM)U= _ 

= u 2 \ x=x O(e N ), (1.6) 

where N > is an arbitrary number. This follows from the fact that the 
difference (Si — S 2 )\ x in parentheses in the exponent is positive. 

Thus, at each point in formula (|1.5p . it is necessary to choose the term 
where the function Sj is minimal. Such a choice leads to an expression of 
the form 

u = e -*( x M e <p(x,t), (1.7) 

where $ = <&(x,i) = mm x {Sj(x,t)}. It is clear that expression (jl.7p is the 
leading term of the approximate solution. 

The corresponding Lagrangian manifold is the following one: 

[P 




X 



Figure 4: 
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The vertical line position is such that the above-mentioned squares are 
equal. 

It is interesting to note that, in the one-dimensional case, there is a direct 
connection between the Hamiltom-Jacobi equation 

St + H(x, S x ) = 

and the conservation law of the form 

d 

u t + —H(x,u) = 0, 

where u = S x and the velocity of the vertical line is defined by Rankine- 
Hugoniot condition corresponding to the conservation law. If the Lagrangian 
manifold has a jump, then the corresponding value function 

$(x, t) = min x S(x, t) 

has a jump in the first derivative. 

Assumption of the actual analyticity for all objects provides that there 
is no "concentration of singularities", — each one can be considered sepa- 
rately in a sense and one can construct a value function as a solution of the 
Hamilton-Jacobi-Bellman equation. 

What about the amplitude function? As was mentioned above, it is a 
solution of the transport equation 

(pot + (x, Vv? ) + ^tr(%(x, VS)S XX = 0, (1.8) 

Generally, the velocity field calculated from the Hamilton system has a 
jump simultaneously with a jump in p = VS (and then in x). Thus, the 
problem (still open!) is to solve the transport equation in a discontinuous 
velocity field. We avoid this problem considering the squared solution of the 
transport equation p = y>jj. Madelung, [8] (about 100 years ago!) observed 
that it satisfies the continuity equation 

p t + (V,ip) + ^trP x? = 0. 

in the smooth case. 

Our case is more complicated: we again have a discontinuity velocity 
field. There are few approaches to this problem solution: the theory of 
measure solutions (F. Murat, P. LeFloh, B. Hayes, T. Zang, Y. Zheng et 
al., [TOl [HJ EH D2J H3J Ej ) , the box approximation (M. Oberguggenberger, 
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M. Nedel'kov, [9]), the weak asymptotics method (V. Danilov, V. Shelkovich, 
[291 USB EH [321 [33]), the generalized characteristics method (V. Danilov, 
D.Mitrovich and V. Danilov, [Ml US [281 [26] ) . The last allows one to con- 
struct a solution to the continuity equation in the case where a singular sup- 
port of the velocity field is a stratified manifold with smooth strata which 
are transversal to the (incoming!) trajectories of the velocity field. 

If for some time interval [t±, t2], the singular support of the velocity field 
preserves its structure (the mapping of the singular support induced by the 
shift along the Hamilton flow is a diffeomorphism) , then it is possible to 
show that the singular support of the velocity has the required structure. If 
the structure is changing (e.g., a jump appears, see Fig.5 and Fig.4-the last 
step of evolution in time) , then one can use the weak asymptotics method to 
construct a global solution to the Hamilton-Jacobi and continuity equations. 
This approach is based on a "new (generalized) characteristics" constructed 
by V. Danilov and D. Mitrovic, [241 \25\ [26] in the case where the strata of 
the singular support are of codimension 1. 




Figure 5: 

The main idea of this approach is to consider the singularity origination 
as a result of nonlinear solitary wave interaction. 

A simple example is the Hamilton flow corresponding to the heat equa- 
tion from the previous example. The Hamilton-Jacobi equation in this case 
is equivalent to the Hopf equation for the momentum p: 

Pt + (p x f = 0. 

The solution in this case has the form 

p = p + a(H((f)i - x){4>i - x) - H(<fi 2 - x)(4>2 - x)) 

and is plotted below. 



S 



Figure 6: Classical characteristics 




Figure 7: New characteristics 
To consider p 2 , we have to calculate the product 

Here the following equality holds: 



H(<f>! - x)H(<f> 2 -x) = BiP—Z±)H{<h. - x) 



+(1 _ B (**—^))H(<h -x) + D ,(fi), 




Po 



> 



() 



X 



Figure 8: Plot of the function p 



where [i is an arbitrary small parameter and Orj'(p) is a small quantity in 
the sense of distributions, 



for each t/j, which is a test function. The time evolution of the function p is 
such that the slanting intercept of a straight line preserves its shape until 
it takes the vertical position and then a jump begins to propagate. This 
means that, at every time instant, the solution anzatz can be presented in 
the form of a linear combination of Heaviside functions. This allows one 
to use a formula which express the product of Heaviside functions as their 
linear combination, and hence, uniformly in time, we see that the functions 
p and p 2 (the last up to a small quantity) belong to the same linear space, 



for detail, see [331 [Ml 113 [26] . 

Thus, we can prove the following theorem. 

Theorem 2. Assume that the following conditions are satisfied for t G [0, T], 



form a stratified manifold with smooth strata and Hess^P(x,^) > 0. 

Then there exists a generalized solution p of the Cauchy problem for 
continuity equation in the sense of the integral identity introduced in \31\ \32\/ 
and at the points where the projection it is bijective, the asymptotic solution 
of the Cauchy problem for Kolmogorov- Feller type equation has the form 



(0 D ,(n),^) = 0( f j.) 



T > 0: 



(1) there exists a smooth solution of the Hamiltonian system, 

(2) the singularities of the velocity field 

u = V f P(x,VS) 
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2 Backward in time motion 



As was shown above, all that we need to go forward in time is the Hamilton 
system: 

x = V^P(x,p), x\t=o = a, (2.1) 
p = -V x P(x,p), p\ t =o=p(a). 
Let us change the time direction as t — > — t, then 

-X = V^P(x,E), (2.2) 

-E = -V x P(x,E),. 
We want to solve the inverse problem: 

X(a,0) = x(a,T), E(a, 0) = p(a, T) 

The right-hand sides are considered as given data, and we are looking for 
X(a, t), E(a, t) for < t < T. Obviously, in our case the solutions have the 
form 

X(a,t) =x(a,T-t), E(a,t) =p(a,T-t). 

Conclusion: we can use the "same" trajectories to move forward and back- 
ward in time. But the incoming trajectories become outcoming and vice 
versa. 

Corollary 1. Stable jumps become unstable. 

But if there are no jumps (singularities of the projection mapping ir : 
A f n — > M™), then our geometry (and the asymptotic solution!) is invertible 
in time. This means that if we take the Cauchy problem solution u for 
parabolic PDE such that 

u\ t =o = exp (-S (x)/e)ip°(x), 

then the asymptotic solution for t = T has the "WKB" form 

uJ\t=T = exp(-S(x,T)/e)tpo(x,T). 

Then taking the last function as the initial data for parabolic PDE in inverse 
time (let v as (x,t) be its asymptotic solution), we get: 

Vas (x,T) = exp {-S (x)/e)<p°(x)(l + 0(e)) 
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It can be easily verified in the case P(x, £) = £ 2 i.e. for simplest heat 
equation. If one constructs the solution of inverse heat equation with initial 
data at t = T of the form 



■Uas|t=T = exp(-S(x,T)/e)ip {x,T) 

using the Green function and calculates the integral be saddle point method 
at t = Othen the following result will be obtained 



n as | i= T * G \t=o = exp (-5 (x)/e)/(x)(l + 0(e)) 



We want to stress once again that this statement is true if there is no sin- 
gularities of projection mapping tt : A l n — > R™. 
But a jump brings problems: 





Figure 9: 



There is no unique reconstruction of the part of Lagrangian manifold 
coming to the vertical line as t increases (see Fig. 9)! But, fortunately, we 
can move ahead using the sense considerations. The main point is that the 
function S cannot attain its minimum (maximum) inside the "terra incog- 
nita" . This allows one to calculate the integrals containing the reconstructed 
solution without taking "terra incognita" into account in the case where the 
integrand support contains this "terra incognita", and we can formulate the 
following statement. 

Theorem 3. Let the symbol P(x,£) defined by (jl.ip be such that the 
function B and measure [x do not depend on x and V = 0. Let u £ (x,t) 
be a solution of the Cauchy problem to a Kolmogorov-Feller-type equa- 
tion, and assume that, for some t G (0,T), there exists a logarithmic limit 
S(x,t) = lim e _s.o ( — elnu e ), namely, the action function and the generalized 
amplitude p = w — lim^o ex.p(2S / e)v% . 
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Let the singular support of S(x,t) be a stratified manifold and p re g be 
the regular part of the generalized solution to the continuity equation in the 
distribution sense. 

Then, for an arbitrary test function <p = (j){x) > from the Schwartz 
space and for all t £ [0, T\, the limit as e — > of the integral 

cj(u £ - pV* eM-S/s))(/>dx (2.3) 

is equal to 0, where C~ l = f u £ (pdx. 

This result can be extended to C^-test functions and to the set of smooth 

1/2 

functions integrable with weight p r i g exp(— S/e) under the same assumption 
as in the theorem above. 

Let Q(t) be a subset in i?™, where we cannot define the functions S and 
(p uniquely by the characteristics. Here we have two possible statements: 

(i) let £l(t) G supp0, then relation (|2.3p is true; 

(ii) assume that x 6 R , Q(t) is a union of segments Ik, and the intersection 
between some 1^ and supp^ is not empty but 1^ does not belong to 
ft(t), then the limit of the integral as e — > 

exp(^/e) J (u e - pl^exp(-S/e))<f>dx 

equals 0, where ^ is the minimal value of S at the ends of Ij,. 

These statements actually mean that, from the viewpoint of the weak 
sense (momentum), the density reconstructed arbitrarily inside the "terra 
incognita" and according to the characteristics outside it can be used in the 
same manner as the leading term of the asymptotic solution constructed 
earlier by Maslov's tunnel canonical operator and its modifications, [U [191 

I22JE31IZI1I2D]. 

Now I will briefly speak about the proofs of the statements about the 
invertibility in time that was formulated above. Each of them can be divided 
into two parts. First, it is to prove that, for all smooth reconstructions 
of Lagrangian manifold in a "terra incognita" domain, the corresponding 
function So = Sq{x) (for some fixed t) cannot attain its minimum value 
inside this domain, see the lemma below. This allows one to apply the 
Laplace method for calculating the integrals mentioned in those statements 
taking into account that, due to this method and Lemma, the results of 
these calculations do not depend on the values of the integrands inside the 
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"terra incognita" domain. The proof is finished by taking account of the 
estimation 

u as -u = 0(e)u as , 

which is true outside the singular support of the function S, see Theorem 
above. Now I formulate the Lemma. 

Lemma 1. Let the symbol P(x,^) defined by (jl.ip be such that the function 
B and the measure /i do not depend on x and V = 0. Then the function S, 
i.e., the action function corresponding to the Lagrangian manifold, cannot 
attain the minimal value inside the "terra incognita" domains. 

We begin the consideration with a particular case when operator symbol 
P(x, £) does not depend on x and restrict ourselves by one dimensional case 
studying. Let (a,b) is an interval inside the "terra incognita" domain, and 
let x~o £ (a, b). We will proceed by contradiction. Assume that the function 
S attain its minimal value at the point xq S (a, b). We prove that, along 
the trajectory of Hamilton system whose projection starts at xq S (o, b), the 
following inequality is true: 

Hi—" 40 - <2 - 4) 

This inequality leads to a contradiction because of the assumption that 
xq belongs to the "terra incognita" domain, and hence it belongs to the 
projection of the image of the singular (vertical) part of the Lagrangian 
manifold under backward in time shift along the trajectories of the Hamilton 
system. In turn, this means that the projections of all trajectories whose 
starting points are projected to the "terra incognita" must intersect at a 
point for the forward in time motion. So the above-mentioned inequality 
leads to a contradiction. To prove this inequality, we write the projection 
of the Hamilton system trajectory starting at xq. It has the form 

x = x +tP^(p ), Po = S 0xo (x ). (2.5) 

It is clear that po = Sq Xo (xo) = 0. Thus, 

dx 

0^-\x =x = 1 + tPtf(Q)S 0xoXo (x ). (2.6) 
Taking into account that that 

d 2 P 

W > (2 ' 7) 
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because of the convexity of P and 



§ >- < 2 - 8 > 

due to the assumption that xo is a point of minimal value, we get the needed 
inequality (|2.4p . A multidimensional case differs from the case considered 
by changing the scalar values P% and Sq xo by vectors (gradients). In turn, 
this gives matrix inequalities in (|2.7p and (|2.8p . The problem is to prove that 
the eigenvalues of the matrix P^(0)So xoxo (xo) are nonnegative by using (12. 7j) 
and (|2.8p . For this, we can make a change of variables reducing the matrix 
So xqxo (xq) to diagonal form. This transformation induces the corresponding 
transformation in the p-plane that transforms the matrix P pp (0) to a new 
symmetric positive matrix. Now we note that the principal minors of the new 
matrix product are products of matrix-factor principal minors (because the 
second one is of diagonal form) . The determinants of the matrix-factor prin- 
cipal minors are nonnegative, so the spectrum of the matrix P^(0)So XQXO (xo) 
is also nonnegative and we again get (|2.4p . To finish our consideration, we 
have to investigate the case where the symbol P = P(x,£) depends on x. 
We again stay at the point xq, where the function Sq attains its minimum. 
If so, then 

d S (x ) , 
dx 

and, by assumptions, p = for t > 0. 

The system for the matrices and follows from the Hamilton 
system and has the form 

d dx d 2 P dx d 2 P dp 

+ ^7^- ( 2 - 10 ) 



p\t=o = —^z l*o=xo = °i ( 2 - 9 ) 



dt dxo d^dx dxo d£? dxo 

d dp d 2 P dp d 2 P dx , n ««s 
£_ = £_ _| (2.11) 

dt dxo dxd£ dxo dx 2 dxo 

Because of our assumptions (see the Lemma formulation), we have P\^=o = 
0, Pxx\z=o = an d P£x\z=o = 0- This means that, along the Hamilton 
system trajectory starting from the point x = xo, p = 0, we have p = and 
jji-j^r = 0. Thus, along the above-mentioned trajectory, equations (|2,10p 
and (|2.1ip have the form 

d dx d 2 P(x ,0) dp 

dtdxo di 2 dxo' [ ' ' 
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d dp 

dt 8xq ^ 

Integrating over t, we can transform the first equation to the form of Eq. (|2.6I) 
and apply all above arguments concerning this equality. 
This note ends the proof. 

The statement of the Lemma can be generalized as follows. Let B be a 
linear function in x, 

n 
k=l 

Then, instead of (|2.13p . we get 

d dp dB dp 



dt 8xq dx dxo 

*1 f\ fTI f~\ 1 t~ 

dx 



where the matrix jM- has diagonal form with elements equal to constants. 



Then 

dp t 8B_ dpo 

- — = exp a* - — 
dx oxq 



and 



dx „ d 2 P /■* t'pt dp 
- E + — - / exp 9 *o dt' • 



dx <9£ 2 J dx 

= J exp*' If dt' y(^J exp*' Weft' 



-1 ^2 



| d 2 Pd Po 
d^ 2 dxo 



Here we used the fact that the matrix 



/"* t'22. , / 
/ exp 1 ex dt 

Jo 



is in diagonal form with positive eigenvalues. 

Thus we came to the relation with the same properties as (|2.6h . One 
can proceed further and generalize the statement to the case of an arbitrary 
drift. But up to now the presence of the potential destroys our picture and 
I will think about it. 
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